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Abstract 

Imposing initial conditions to nonequilibrium systems at some time to leads, in renormal- 
ized quantum field theory, to the appearance of singularities in the variable t— to i n relevant 
physical quantities, such as energy density and pressure. These "initial singularities" can 
be traced back to the choice of initial state. We construct here, by a Bogoliubov transfor- 
mation, initial states such that these singularities are eliminated. While the construction 
is not unique it can be considered a minimal way of taking into account the nonequilibrium 
evolution of the system prior to to- 
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1 Introduction 



Nonequilibrium dynamics in quantum field theory has become, during the last years, a 
very active field of research in particle physics [1-6], in cosmology [7-20] , and in solid state 
physics plR . The outline of typical computational experiments is as follows: a quantum 
field i/)(xi,t) is driven by a classical field degree of freedom (Higgs, inflaton, condensate) 
<fi(t) which takes an initial value away from a local or global minimum of the classical or 
effective action; the time development is then studied including the back reaction of the 
quantum field in one-loop, Hartree or large-N approximations. The initial state of the 
quantum field ip is usually taken to be the vacuum state corresponding to a free field of 
some "initial mass" m(t ) or a thermal state built on such a vacuum state. In Friedmann- 



Robert son- Walker (FRW) cosmology |TT| , 18, 19 1 the usual initial quantum state is chosen 



to be the conformal vacuum, again corresponding to the initial mass m(to). While such 
choices seem very natural, they are not necessarily appropriate; nevertheless, this point 
has received little attention up to now. The reason why we address this question is the 
occurrence, in some dynamically relevant physical quantities, of singularities in the time 
variable which are related to the choice of initial state. We will in fact show that these 
singularities can be removed by more appropriate choices. 

The origin of these singularities can be traced back to a discontinuous switching on of 
the interaction with the external field 4>(t). This interaction is given by a time-dependent 
mass term [] m 2 (t) = m 2 (0) + V(t) where for the simplest case of a A$ 4 theory H [EJ 

In FRW cosmology the (conformal) time dependent mass term reads 0, [19 



MV) = ^ + ( f -i) BW + gl}}. (1.2) 

In this case M 2 (r) and therefore V(r) = M 2 (r) — M 2 (0) also contain the scale parameter 
a{r) and the curvature scalar R(t). 

A free field theory vacuum state corresponding to the mass m(0) would be an appro- 
priate equilibrium state if V(t) stayed zero for all times. However, at t = the potential 
V(t) changes in a nonanalytic way. This is unavoidable since at least the second derivative 
of (j)(t) becomes nonzero on account of the equation of motion. As a result of these discon- 
tinuities relevant physical quantities develop singularities in the time variable t at t — 0. 
In the case of A$ 4 theory in Minkowski space such singularities only occur in the pressure. 
In FRW cosmology the problem becomes more acute. There, even the first derivative of 
V(t) necessarily becomes nonzero at t = 0; indeed, even with a constant external field 
the initial state could not be at equilibrium; this manifests itself by a nonvanishing first 
derivative of the scale parameter induced by the Friedmann equations. Furthermore, in this 
case both energy and pressure become singular; since they enter the Friedmann equations 
this singular behavior also affects the dynamics. 

4 We choose to = for convenience. 
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Singularities arising from imposing initial conditions to quantum systems have been 
noted for the first time by Stiickelberg [p3] , who called them "surface singularities" ; they 



are briefly mentioned in the textbook of Bogoliubov and Shirkov pi |. The "Casimir 



effect" arising from initial conditions has been discussed by Symanzik |25| . In the context 
of nonequilibrium dynamics in FRW cosmology the occurrence of such singularities has 
been noted by Ringwald 0. In the following we will refer to these singularities as "initial 
singularities" . 

Imposing an initial condition at some time to does not mean, in most applications, that 
one assumes the system to have come into being at just this time. Rather, t is usually 
chosen as a point in time at which one can make, on some physical grounds, plausible 
assumptions about the state of the system. Clearly, if the system is not at equilibrium 
after to, it will not have been so before. Therefore, the initial state should take into account, 
at least in some minimal way, the previous nonequilibrium evolution of the system. Such 
a minimal requirement is the vanishing of initial singularities. It is the aim of this paper 
to specify such initial states. 

The plan of the paper is as follows: in section 2 we present the basic equations and 
formulate the problem for the case of a A$ 4 theory; in section 3 we discuss an appropriate 
choice of the initial quantum state such that the singular behavior in the time variable is 
removed; the modified renormalized equations for the nonequilibrium system are given in 
section 4; we end with some concluding remarks in section 5. 

2 Formulation of the problem 

We consider a scalar A$ 4 theory without spontaneous symmetry breaking. The Lagrangian 
density is given by 

£ = -d.Qd^ - -m 2 $ 2 - -$ 4 . (2.1) 
2 m 2 4! y J 

We split the field $ into its expectation value and the quantum fluctuations ip' 

$(x,t) = 0(t)+^(x,t), (2.2) 

with 

0(t) = ($(x,t)) = ^^, (2.3) 
where p(t) is the density matrix of the system which satisfies the Liouville equation 

i^- = [n(t),p(t)}. (2.4) 
The Lagrangian then takes the form 

C = Co + & , (2.5) 
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with 



4! 6 4 6 



The equation of motion for the field <j)(t) is given by || 

4>(t) +m 2 0(t) + ^0 3 (t) + -^0(t)G ++ (t,x;t,x) = . 
6 z 2 

Here G ++ is the ++ matrix element of the exact nonequilibrium Green function 
in the background field 4>(t). For a pure initial state \i) it can be written as 

-zG ++ (t,x;t',x') = (i\Tip(t,x)ip(t',x')\i) . 



(2.6) 
(2.7) 

(2.8) 



(2.9) 



If the classical field is spatially uniform the equation of motion for the field ip(t, x) is given 
by 



§-A + - 2 + ^ 2 (t) 



V(t,x) = 



We introduce the notations 



m 2 (t) 
u k (t) 



and 



m 2 + ^<j) 2 (t) 
k 2 + m 2 (t) 



k 2 + mt 



We will discuss the choice of mo below. We define the 'potential' V(t) as 

V(t)=u 2 (t)-u 2 

We further introduce the mode functions for fixed momentum Uk{t) exp(zk ■ x) 
satify the evolution equation 



at 2 



U k (t) = ; 



we choose the initial conditions 

U k (0) = 1 U k (0) 



-lUJ k0 



(2.10) 

(2.11) 
(2.12) 

(2.13) 

(2.14) 
which 

(2.15) 
(2.16) 



J Note that the functions U k (t) depend only on the absolute value of k. 
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The field ip can now be expanded as 

^(t,x) = J {2n p 2uko [a(k)U k (t)e^ + J(k)U* k (t)e-^} , (2.17) 

where the operators o(k) satisfy 

[a(k),at(k')] = (27r) 3 2^ 3 (k-k') (2.18) 

If the initial state \i) is chosen as the vacuum state corresponding to the operators a(k), 
i.e., as satisfying a(k)|i) = 0, we obtain the Green function G ++ (t, t'; x — x') as 

d 3 k r „ ,_sr^,.,s„,. „ / .is , .s „/ ., .s, 3 jk(x-x') 



(2.19) 

The Green function at equal space and time points then reads 



The resulting equation of motion for the classical field <f>(t) is 

4>(t) + m 2 0(t) + ^0 3 (t) + ^(f).F(f) = 0, (2.21) 
where we have introduced fluctuation integral 

^=i<S^ m)r (2 - 22) 

It determines the back reaction of the fluctuations onto the classical field <f>(t). 

We further consider the energy density and the pressure. The energy density is given 

by 

£ = \^) + V(m + . (2.23) 

Calculating the trace over the Hamiltonian for the same initial state we obtain 

Using the equations of motion it is easy to see that the time derivative of the energy density 
vanishes. 

The pressure is given by 

' = ^ + 1 p|k K + (()|2 + t^A - £ ■ < 225 > 



While the pressure does not enter the dynamics here, so it does in FRW cosmology. 

The expressions for the fluctuation integral, the energy density and the pressure are 
divergent and one has to discuss the renormalization of this theory. We have presented 
recently p2| a fully renormalized framework for nonequilibrium dynamics. The main tech- 
nical ingredient of this analysis is the perturbative expansion of the functions U k (t) with 
respect to orders in the potential V(t). We write the functions U k as 



U k {t) 



-iuj ko t 



1 + h k {t)\ 



and expand further in orders of the potential V(t) as 



n=l 



We also introduce the partial sums 



fc£°(*) = E^ ; (*) 



(0, 



l=n 



(2.26) 
(2.27) 

(2.28) 
(2.29) 



so that 

h h {t) = h£\t) = hP + h { p(t). 

The integral equation for the function h k (t) can be derived in a straightforward way from 
the differential equation satisfied by the functions U k (t); it reads 



h k (t) 



2uko Jo 



We obtain 



h 



(i) 



dt' (e^ ^ - l) V(t') [1 + h k (t')} . 



dt' (e 2 ^"(*-*') - l) V(t') . 



2u> k0 Jo 



(2.30) 



(2.31) 



Using integrations by parts this function can be analyzed with respect to orders in u k0 via 



h ( k \t) 



2u; kQ Jo 



dt'v(t')+j2 

1=0 



fcO 



1+2 



V®(t) -e 2!UfcO V (!) (0) 



-i \ n + 2 rt 



(—) 



dt'e^^V^it') 



[2.32) 



where V®(t) denotes the Ith. derivative of V(t). For energy density and pressure we need 
the expansion of h k l \t) as well. From Eq. Q2.32D and the relation 



2iuj ko h k 1) 



dt'Vit') 



(2.33) 



we find 



h ( k \t) 



E 

1=0 



( 



\2u 



k o 

, \™+! rt 



2u. 



fcO 



V®(t) -e 2iu > kot V (l) (0)] 

dt f e 2iu; k0 (t-t') V (n+l)^ 







(2.34) 



In the following we will need the real and imaginary parts of this expression; we introduce 
the following useful notation: 

C(f,t) = f dt'f{t') C os{2u m t') , (2.35) 
J o 

S(f,t) = f dt'f{t')sm{2u m t') . (2.36) 
J o 

We now insert the perturbative expansion into the fluctuation integral to obtain 

/d ^ k ( 

d 3 k f V(t) V(0) tn , V(0) . tn , V(t) 
+ TtV cos(2a; fc ot) + ~V sm(2w M t) + 



(2vr)32a; fc o [ 2^ v KU ' 4^ v - ' 



fcO 



1 cos(2^ot) - ^ r C(V ; ,t) + 2Re4 2) + l^fcl 2 ! • (2.37) 



The first two terms in the parenthesis of the second expression, i.e., 1 and V(t)/2u>l lead 
to divergent integrals which have to be absorbed by the renormalization procedure. This 
has been discussed in [52]. There, the mass m was chosen to be the 'initial' mass m(0) 



r d 3 k V(0) 1 , . , 

/ TTTT^o TtV cos 2a; fc0 * - -t-j In 2m t as t -> . 2.38 



(see (|2.11 )). It was shown that the renormalization counter terms do not depend on this 
mass but can be chosen to contain only the renormalized mass m corresponding to the 
perturbative ground state at <fi = 0. With this choice of initial mass V(0) is zero and the 
fluctuation integral is nonsingular at t — 0. 

If, on the other hand, we choose mo = m it is obvious that the divergencies are absorbed 
by the counter terms depending only on the perturbative mass m, but we are faced with 
an initial singularity arising from the third term via (see Appendix B) 

d 3 k V(0) tn , 1 

fo^ao ^~2~ cos ( 2uJ kot) ^ - — 

(2'7i)- i 2uj k0 2uif, 8n z 

Of course, nobody has made such an 'unnatural' choice of the initial mass, this initial 
singularity can be avoided trivially by choosing m = m(0). It is important to note, 
however, that the renormalization can be performed in a way independent of the initial 
condition in both cases. The difference between the two approaches is in the initial 'vacuum' 
state. These different initial states are related by a Bogoliubov transformation (see also 
Appendix A). So Bogoliubov transformations can be used to avoid initial singularities. 

In (|2.37|) we have extended the expansion of Re/4 (0 to display also the terms of 
order uj^ which depend on V(0) and V(0). These terms do not lead to divergencies in the 
fluctuation integral; however, in the energy and pressure they appear multiplied with u>l 
and/or k 2 . While the energy stays finite the pressure behaves in a singular way via 



Pfluct. 



smg 



d 3 k 



(2vr) 3 2^ 



\4ux, 



sm(2u k0 t) - —\- cos(2uj k0 t) \ . (2.39) 



J k0 OU7 fcO ) 

The behavior of the momentum integrals is given in Appendix B, they result in a l/t 
singularity proportional to V(0) and a logarithmic one proportional to V(0). Therefore, 
these terms have to be removed as well. 
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3 Removing the initial singularity 



We have seen in the previous section that nonzero initial values of V{t) and its derivatives 
lead to initial singularities. The clue for dealing with these terms has already been indi- 
cated: the leading singularity can be removed by a Bogoliubov transformation from the 
perturbative vacuum to a vacuum corresponding to free quanta of the initial mass m(0). 
We expect, therefore, that the other singular terms can be removed in this way as well. 
We define a general initial state by requiring that 



a(k) - p fc a f (k) 







(3.1) 



The Bogoliubov transformation to this state is given in Appendix A. If the fluctuation 
integral, the energy and the pressure are computed by taking the trace with respect to this 
state the functions Uk(t) are just replaced by 



F k (t) = cosh( 7fc )t4(t) + e Lhk smh^ k )U* k (t) , 
where 7^ and 5 k are defined by the relation 

pk = e tS tanh(7 fc ) . 
The fluctuation integral now becomes 

d 3 k 



(3.2) 



(3.3) 



(2tt) 3 2c^ 
d 3 k 



(27r) 3 2^ 
Expanding as before we find 
d 3 k 



Fk{t)? 

{cosh(2 7fc (t))|f4(t)| 2 + sinh(2 7fc )Re (e^U^t))} 



(3.4) 



(2tt) 3 2^ l 
V(t) V(0) 



cosh(27^ 



1 - 



V(t) , V(0) 



2^o 



8w fc0 



cos(2ukot) 



cos(2u;a;o£) + 



-C(V,t) + 2Re4 2) + \h k 



V(0) 



sin(2a;fcot) 



feO 



k0 



+ sinh(27 fe ) cos(S k ) cos(2c^. t) 

- sinh(2 7fc ) sin(4) sm(2uj k0 t) + sinh(27 fc )Ree- 2 ^ fc0 '~ i<5 (2h k + H 



(3.5) 



Let us first discuss how to get rid of the most singular term, proportional to V(0). Requiring 
this term to be compensated by the terms proportional to sinh(27fc) we find 



tanh(27 fc ) 



= 0, 



^(0) 
2^o 



(3.6) 
(3.7) 
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As explained in Appendix A the standard Bogoliubov transformation from the perturbative 
vacuum with mass m to the vacuum corresponding to quanta with the initial mass m(0) 
is mediated by a function j' k (k) satisfying 

which implies 

tM) = 2^(0) • < M > 

We see that 7^ and j' k agree asymptotically to leading order in l/u k0 . So requiring that 
the most pronounced initial singularity vanishes leads essentially to the usual choice for 
the initial state, namely mo = m(0) and therefore ^(0) = . The analysis of subleading 
terms in the difference between 7^ and 7^ becomes somewhat cumbersome. After we have 
convinced ourselves that the Bogoliubov transform is the right technique for getting rid 
of initial singularities we will therefore choose m = m(0) as everybody does and apply 
this technique to get rid of the remaining singularities. So from now on 1^(0) = and 
u k o = (k 2 + m 2 (0)) 1//2 . Requiring that the terms proportional to V^(0) and V^O) vanish 
leads to the conditions 

tan(5 fc ) = 2 ^o^y, (3.10) 



tanh(27, 



k 



t2i 



V 2 (0) + 



V 2 (0) 



J k0 



(3.11) 



Using these functions we are now ready to formulate the renormalized equation of motion 
and the energy momentum tensor. 



4 The renormalized equations 

We have given the bare equation of motion and energy momentum tensor in section 2. The 



renormalization for the original initial state has been discussed in ||22|| . We have to ensure 
now that the scheme used there is not spoiled by the improved initial state. The main new 
feature in the fluctuation integral, the energy density and the pressure is the appearance of 



factors cosh(27fc) and sinh(27fc). We will need their asymptotic behavior. Using Eq. ( 3.11 
we have 

(4-12) 

The factor cosh(27fc) is equal to 1 for 7^ = 0; we will need the difference 



cosh(2 7fc ) -1 = 2 sinh 2 ( 7fc ) ^ . (4.13) 

32a; 



AO 
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There are new terms proportional to sinh(27^.); this factor behaves as 

k-,00 \V(0)\ 



sinh(27 fc ) 



4^o 



(4.14) 



The dimensionally regularized fluctuation integral (|3.5| ) takes, after cancellation of the 
singular integrals induced by Eqs. (|3.10| ) and ( |3.1ip , the form 



J- Teg (t) 



"0 (T 4- n v ' t 



d 3 k 



V(t) 
sinh 2 (7 fc ) 



Mo 



cosh(27 fc ) 



V(t) 1 



, , C(V : ,t)+2Re4 2) + |/i fc l 2 



+ sinh(27 fc )Ree- 2 ^ fe0 ^ i ' 5 (2/i fc + h 



(4.15) 
(4.16) 



Here we have introduced the abbreviation 



T 2 4-1 ^ 

L = - + In -— t, 7 



7?); 



(4.17) 







Using the estimates (|4.13|) and ( 4.14 ), and using the fact that the mode function h\. behaves 
as u^q we see that the momentum integral is convergent. 
Introducing the counter term Lagrangian 



C c ,. = Um^ + 



(4.18) 



the fluctuation integral gets replaced, in the equation of motion ( [2.2 1|) for 4>(t), by 

25m 2 S\ 



3~ Em 3~ reg 



With the standard choice 



5m 2 
8X 
L 



[L + l) 



Am 
32tt 2 

32vr 2 ' 

2 , 4vr/i 2 

- + In — — - 7 

e m z 



(4.19) 

(4.20) 
(4.21) 
(4.22) 



jF fin is indeed finite. 
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The calculation of the energy density proceeds in an analogous way. The fluctuation 
energy becomes, using dimensional regularization, 



= ^/(2^k^ i(t)|2+( ^ +l/(t))|Fi(t)|2 } 



t0 (L + 3 ) V{t) (L+l\ V2{t) L 



+ 



1 r d 3 k 

2 J (2txY2u 



2 J (2tt)32^ 
+V(t) cosh(2 7fc ) 



2sinh 2 ( 7fc ) 



2^o + V(t) 1 



+V(t) sinh(27 fe )Ree- 2 ^ fcl '^ i,5fc (2h k + /i 
+ cosh(27 fe )|/i fc | 2 



+ sinh(2 7fc )Ree- 2 ^- i5fc \h 2 k - 2iu k0 (l + h k )h k 



(4.23) 



(4.24) 



The divergent parts are cancelled by the counter terms 

Set. = 5A + ^m 2 2 (t) + ^0 4 (t) 
with the 'cosmological constant' counter term 

6A = ^ L + J- 



(4.25) 



(4.26) 



Finally we have to consider the pressure. We find for the regularized fluctuation part: 



Pfluct = — £| 



fluct 



< v(t) 2 i>(t) 



._. 777, — 

96tt 2 48tt 2 96tt 2 



+ 



/ 



(2ir) 3 2u k0 



sinh 2 (7 fe ) 



k 2 



V(t) , V(t) 



+ 



+ cosh(27 fe ) 



c(v-',t) 

&4> 



+ 2Re4 2) + 



+ cosh(27 fc )|/i fc | 



(4.27) 



sinh(27 fc )Ree- 2 ^°*- i,5fc 



A: 



- ^lo + "o") ( 2/lfe + ^) + ^ ~ 2ttj fc0 (l + 



In order to cancel the divergent term proportional to V one introduces a counter term for 
the energy momentum tensor 



ST, U = A(g^d a d a - d^)& , 



which leads to a counter term 



cP 2 A 

Pc , t= A^( t) = -V(t). 



(4.28) 



(4.29) 
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in the pressure. We choose 

* = -«£i £ - < 4 - 3 °> 

The remaining momentum integral is finite and nonsingular in t. For most of the terms 
this can be seen by inspection using Eqs. ( |4.13| ), ( |4.14j ), and the expansions ( |2.32| ) and 
( |2.34j ) . There are some internal cancellations which are, however, the same as for the case 
7fc = already discussed in [E2|. The only new, potentially singular term is 



The leading singular behaviour is given by 

V(0) / dt'V{t') / ; -V( - cvl + -) cos(2^ t) . (4.32) 

While the momentum integral behaves as In t as t — > the time integral behaves as t 2 since 
V(0) = 0. So the renormalized pressure is indeed nonsingular at t — 0. 

From the analysis of divergent integrals given in this section it is obvious that only 
the leading asymptotic behavior of 7& is relevant, more precisely, only the terms of order 
u>kQ and u>k Q 4 . This means that any Bogoliubov transformation whose function ■jk has this 
leading asymptotic behavior is equally suitable for defining an appropriate initial state. 

We have formulated our modified renormalized equations for A$ 4 theory in flat space. 
The generalization to a scalar field in a flat Friedmann- Robert son- Walker universe is 
straightforward and the cancellation of singular terms in the energy density and the trace 
of the energy momentum tensor proceeds in the same way. 



5 Conclusions 

We have considered here the choice of initial states for a nonequilibrium system in quan- 
tum field theory. Our considerations arose from the problem that logarithmic and linear 
singularities in the variable t — t Q appear in the energy momentum tensor and affect the 
dynamics of FRW cosmology. We consider such singularities - and their consequences - as 
unphysical, at least if t is just some conveniently chosen point in time within a continuous 
evolution of the system. Most authors, including ourselves, have chosen initial states that 
correspond to equilibrium states of the system. We have constructed here improved initial 
states for nonequilibrium systems in such a way that the appearance of initial singularities 
is avoided. These states are obtained from the usual 'vacuum' states by a Bogoliubov 
transformation. The essential part of this transformation is a Bogoliubov 'rotation' of the 
creation and annihilation operators at large momentum. The construction presented here 
specifies a transformation of creation and annihilation operators at all momenta. It is not 
unique in the sense that it may be arbitrarily modified at small momenta. This non unique- 
ness is, however, nothing else as the freedom for choosing an 'arbitrary', pure or mixed, 
initial state. Our construction can be considered as formulating a minimal requirement for 
choosing such states in the sense that it specifies the initial state of the high momentum 
quantum modes. 
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A Bogoliubov transformation 

In this Appendix we briefly recall some basic features of the Bogoliubov transformation 
(see , e.g., |[28|| ). We start with a vacuum state defined by 



a(k)|0) = . 



(A.l) 



We would like to obtain a new state |0) that is annihilated by a(k) + p^a^k), where pk is 
some complex function of k, i. e. we require 

a(k)-p fc at(k)]|0)=0. (A.2) 

Such a state can be obtained from |0) by a Bogoliubov transformation 

|6>=exp(Q)|0>. (A.3) 

Using the general relations given in |28[] one finds the explicit form of the operator Q as 



^ 2 



d 3 k 



■Ik 



e iSk a){\L)a\-\i) - e^ fc a(k)a(-k) 



(2tt) 3 2^ 

Here jk and 5k are defined by the relation 

pk = e J(5fc tanh7 fc 



so that (|A.2| ) can also be written as 

a(k)|0) = [cosh( 7fc )a(k) + e i5k sinh( 7fc )a t (-k)l |0) = . 



(A.4) 



(A.5) 



(A.6) 



A special class of new 'vacuum' states |0) is obtained when the new creation and annihila- 
tion operators refer to free particles with a different mass ifiQ. In the field expansion ( [2.17| ) 
this means that the energy Uko = (k 2 + m 2 ,) 1 / 2 is replaced by Uko = (k 2 + m 2 ,) 1 / 2 . In this 
case 



l^kO 



5(k) = W^a(k) + J^a+(k) 

V ^fcO V U k0 



and therefore 



Ik = - In — 

L UkO 



while 5k = 0. For k ^> m ,rh the function j k behaves as 

Ik 



1 ~ 2 

77T.Q — 77T.Q 



4fc 2 



(A.7) 
(A.8) 

(A.9) 
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B Some singular integrals 

The singular behaviour in time arises from the following integrals 

r d 3 k 1 

= / To \3o 2-cos(cj fe0 t) (B.l) 

and 

h(t) = I f t , -f - sm(^ot) • (B.2) 

The first integral can be rewritten as 

-r / \ If 00 duj^uJ 2 - m 2 . , 
h W = 4^1 — 7 C0S(2W() 

1 f 00 / 1 771 2 \ 

= / <M / 2 2 COS ( 2 ^) - 2/2 2 C ° S ( 2 ^) • ( B - 3 ) 

The integral over the second term is nonsingular; the first term yields a Bessel function 
Y (2mt), explicitly 

hit) = -^Y (2mt)+O(t 2 ) ~° --i-ln(2mt) . (B.4) 



The integral ^(i) is simply given by 

h(t) = (B.5) 

and therefore 

(B.6) 
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